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ABSTRACT 

The evolution of the temperature in the intergalactic medium is related to the reion- 
ization of hydrogen and helium, and has important consequences for our understanding 
of the Lya forest and of galaxy formation in gravitational models of large-scale struc- 
ture. We measure the temperature-density relation of intergalactic gas from Lya forest 
observations of eight quasar spectra with high resolution and signal-to-noise ratio, using 
a new line fitting technique to obtain a lower cutoff of the distribution of line widths 
from which the temperature is derived. We carefully test the accuracy of this technique 
to recover the gas temperature with a hydrodynamic simulation. The temperature 
at redshift z =(3.9, 3.0, 2.4) is best determined at densities slightly above the mean: 

=(20200 ± 2700, 20200 ± 1300, 22600 ± 1900)K (statistical error bars) for gas den- 
sity (in units of the mean density) =(1.42 ± 0.08, 1.37 ± 0.11, 1.66 ± 0.11). The 
power-law index of the temperature-density relation, defined by T = T^(Ag/A^)'''~"'^, 
is 7 - 1 =(0.43 ± 0.45, 0.29 ± 0.30, 0.52 ± 0.14) for the same three redshifts. The 
temperature at the fixed over-density A = 1.4 is ri.4=(20100 ± 2800, 20300 ± 1400, 
20700 lb 1900)K. These temperatures are higher than expected for photoionized gas in 
ionization equilibrium with a cosmic background, and can be explained by a gradual 
additional heating due to on-going He ii reionization. The measurement of the temper- 
ature reduces one source of uncertainty in the lower limit to the baryon density implied 
by the observed mean fiux decrement. We find that the temperature cannot be reliably 
measured for under-dense gas, because the velocities due to expansion always dominate 
the widths of the corresponding weak lines. 
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Subject headings: cosmology: observations — intergalactic medium — quasars: absorp- 
tion lines 

1. INTRODUCTION 

The Lyman-a forest absorption in the spectra of quasars provides a wealth of information 
about the properties of the intergalactic medium (hereafter, IGM). There has recently been a lot 
of interest in using the distribution of Doppler parameters of fitted absorption lines, measuring the 
total velocity dispersion of the gas, to constrain the temperature of the IGM (Schaye et al. 1999; 
Ricotti, Gnedin, Sz Shull 2000; Bryan &: Machacek 2000). In this paper we develop a method to 
identify and to fit absorption lines, and to obtain the gas temperature from the distribution of 
the Doppler parameters of the lines. Our algorithm is intended to be simple to implement and be 
applied identically on simulations and observations. 

The temperature of the IGM as a function of density is primarily determined by the balance 
between adiabatic cooling and photoionization heating, once ionization equilibrium with the back- 
ground radiation has been established. However, during the epoch of reionization, the heating 
rate is higher because every atom needs to be ionized once (and the ionization can occur on a 
short time-scale compared to the recombination rate) , and the high opacity of the low-density IGM 
implies that high-frequency photons arc absorbed, delivering a much greater amount of heat for 
each ionization (e.g., Miralda-Escudc & Rccs 1994; Hui & Gnedin 1997; Hachnclt & Steinmetz 
1998; Abel & Haehnelt 1999; Gnedin 2000). Other sources of heating may also contribute, such as 
Compton heating by the X-ray background (Madau & Efstathiou 1999), or photoelectric heating 
by dust grains (Nath, Sethi, & Shchekinov 1999). Constraining these sources of heating is one of 
the two primary reasons why we are interested in measuring the temperature. The other reason 
is the need to make accurate predictions for the statistics of the Lya forest flux in order to con- 
strain cosmological parameters (e.g., Rauch ct al. 1997; Weinberg et al. 1997; Croft et al. 1998; 
Hui 1999; McDonald & Miralda-Escude 1999; Hui, Stcbbins, & Buries 1999; Croft et al. 1999b; 
Weinberg et al. 1999; Croft, Hu, & Dave 1999a; Nusser & Haehnelt 2000; McDonald et al. 2000). 
The temperature-density relation affects the predicted relationship between the power spectrum of 
the transmitted flux and the power spectrum of the initial mass density perturbations (Nusser & 
Haehnelt 2000), as well as the predicted mean transmitted flux (Rauch et al. 1997; McDonald et 
al. 2000), which can be used to constrain the baryon density of the universe. 

Recent Lya forest simulations have shown that, when the structure of the absorption systems 
is adequately resolved, the predicted absorption line widths are smaller than observed if the tem- 
perature of the IGM is determined from photoionization equilibrium, well after reionization has 
ended (Theuns et al. 1998, 1999a; Bryan et al. 1999). To solve the discrepancy the temperature 
apparently needs to be higher. Several authors have presented measurements of the IGM tem- 
perature using different methods, generally flnding values moderately higher than expected from 
photoionization equilibrium (Theuns et al. 1999b; Ricotti et al. 2000; Bryan & Machacek 2000; 
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Schaye et al. 2000). 

Our aim in this paper is to provide a new unambiguous measurement of the temperature, 
making a more exhaustive analysis than in previous work of the model uncertainties that result 
from comparing the observational results with a simulation. We develop a new line-fitting method 
as an alternative to the standard Voigt-profile fitting with line deblending, which is much faster, 
unambiguous, and easy to implement. Our method works by essentially assigning one line to each 
sufBciently deep minimum in the transmitted flux, and measuring the line width and central optical 
depth for each line. The gas temperature at each density is then derived from the distribution of 
line widths at each central optical depth. The systematic uncertainties and model dependence of 
the method used to derive the temperature are carefully analyzed, in a more extensive way than it 
was done in previous work. The new method is applied to observational data and to a simulation 
in exactly the same way, computing error bars due to the variance in our observed sample. 

The main idea of the method to measure the temperature of the IGM was suggested by Schaye 
et al. (1999), Ricotti et al. (2000), and Bryan & Machacek (2000). The probability distribution 
of Doppler parameters, P{B), is characterized by a lower cutoff. Be, where P{B) rises sharply, 
with very few lines having narrower Doppler parameters than this cutoff. The idea is that this 
cutoff is a measure of the gas temperature. In general, absorption lines have both a thermal 
and a hydrodynamical contribution to their breadth; however, for any set of lines with similar 
gas temperature, the narrowest ones will be those where the velocity field along the line of sight 
through the absorber is close to a caustic, so that the variation in the fluid velocity is minimized 
and thermal broadening dominates the observed line width. In fact, it was found by Theuns et al. 
(1999b) that the narrowest absorption lines are primarily thermally broadened. 

A tight relationship between density and temperature in the IGM for gas at low densities, 
where shock heating is not very important, is expected theoretically and is found in numerical 
simulations (Hui &; Gnedin 1997; Theuns et al. 1998). This implies that the narrowest lines at 

a given gas density (corresponding approximately to the optical depth at the line center) are not 
selected to have low gas temperature, but low fluid velocity dispersion. This justifies estimating 
the temperature from the lower cutoff of the Doppler parameter distribution. 

In §2 we briefly describe the observational data and the simulation that we use. In §3 we 
describe our line fitting algorithm. In §4 we demonstrate how the line fitter works by running it on 
spectra from the numerical simulation. In §5 we describe our method for estimating the temperature 
from Doppler parameter distribution, testing the conditions under which the temperature can be 
recovered in a model-independent way. In §6 we use the line fitter on the observational data and give 
results for the measured temperatures. The results are discussed in §7. The Appendix describes 
further details of our line-fitting method. 
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2. THE OBSERVATIONAL DATA AND THE SIMULATION 

2.1. Observations 

We use the same set of eight quasar spectra as in Paper I . These spectra have sufficiently high 
resolution and signal-to-noise ratio to measure the shape of each absorption feature. The pixel noise 
is typically less than 5% of the continuum flux level, and frequently as low as 1%. The velocity 
resolution is 6.6 kms^^ (FWHM) and the spectra are binned in 0.04 A pixels. More details and 
statistics of this data set are given in Paper I and references therein. 

The seven quasars from the Ranch et al. (1997) data set have previously constructed lists of 
regions that are suspected of containing metal lines. Our main results include these regions in the 
spectra because they are not positively identified as containing metal lines. 

In Paper I we defined three redshift bins: 3.39 < z < 4.43, 2.67 < z < 3.39, and 2.09 < z < 
2.67, with mean redshifts z = 3.9, z = 3.0, and z = 2.4. We use these same three bins in this paper, 
which contain approximately the same amount of data. 

2.2. Simulation 

We test the profile fitting code and the procedure to measure the temperature on the output 
of the Eulerian hydrodynamical simulation described in Miralda-Escude et al. (1996) (referred to 
as LIO in that paper). The cosmological model used has Qq = 0.4, $7^ = 0.6, h = 0.65, o-g = 0.79, 
and large scale primordial power spectrum slope n = 0.95. The box size of the simulation is 10h~^ 
Mpc, and it contains 288^ cells. We use outputs from the simulation at z =4, 3, and 2. 

2.2.1. Generation of Simulated Spectra 

Lya spectra are computed for a large number of lines of sight along the box axes. There is 
one free parameter that we can vary when computing the spectra, the normalization of the optical 
depth, which we adjust to reproduce the mean transmitted flux of the observations that we are 
comparing to (the values of the mean transmitted flux are taken from Paper I ). Renormalizing the 
optical depth is equivalent to modifying the intensity of the ionizing background, as long as the 
effect of collisional ionization and the change in the gas temperature caused by the different heating 
rate can be neglected (see Theuns et al. 1998for a test that these effects are in fact negligible). 
The optical depth is then mapped to transmitted flux using F = exp(— r). 

For each line of sight through the simulation (parallel to one of the three axes), we estimate 
the effects of continuum fitting by defining the maximum transmitted flux along the line of sight to 
be the continuum fiux, Fc, and dividing the flux in all other pixels in the line by Fc. We map the 
288 cells along a line of sight onto smaller cells, with their size chosen to match the observations 
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that we want to compare to. Finally we convolve the spectra with the instrumental resolution of 
6.6kms~^, and we add Gaussian noise to each cell with a flux dependent dispersion n{F), which 
is taken from Paper I . 

2.2.2. The T — Ag Relation in the Simulation 

Before we describe the method we shall use to measure the temperature, it will be useful to 
examine the temperature-density relation in the simulation. In this paper we parameterize the 
mean temperature-density relation (hereafter referred to as the T — Ag relation ) as a power-law, 
T = ToAg~ , where the gas over-density is Ag = pg/pg, for the purpose of measuring this relation 
from the data. Although the T — Ag relation naturally approaches a power-law form with 7 — 1 — 0.6 
when the thermal evolution is determined by photoionization heating and adiabatic expansion alone 
(Hui & Gnedin 1997), in general it deviates significantly from a power-law. Figure 1 shows scatter 
plots of T versus Ag for the three redshift outputs of the simulation. The solid lines are power-law 
fits to the range 1 < < 2, to show that the T — Ag relation in the simulation is only roughly 
consistent with a power law. There is a substantial dispersion of the temperature at a given density, 
and the relation between the mean temperature and the density deviates from a power-law. For 
example, using the z = 3 simulation output, the best power-law fits in the restricted ranges of 
density =(0-1, 1-2, 2-3) yield 7 — 1 =(0.15, 0.30, 0.39), with a mean fractional temperature 
deviation around the fits equal to (4%, 11%, 18%). The dispersion is due to shock-heating and 
to the variable expansion or contraction histories of the gas at a fixed density. The reionization 
of Hen occurs near 2; ~ 3 in this simulation, heating the low-density gas to a temperature that is 
nearly constant with density. The parameters of the power-law fits to all three simulation outputs, 
in the range A^ =(1-2), are given in Table 1. 

Because of this deviation from a power-law form of the T — Ag relation , a measurement of 
the temperature Tq and power-law index 7 should be understood only as an approximation to the 
true mean T — A^ relation , near the effective density at which the measurement is made, and this 
effective density needs to specified. We use the following form to present our results in §6: 

T = r,(A,/A,)^-^ , (1) 

where A^ is chosen so that the error bars on and 7 — 1 are uncorrelated. 

2.2.3. Definition of Temperature and Density at Points in Spectra 

Each pixel in a spectrum receives optical depth contributions from an extended stretch of 
real space along the line of sight, so no unique temperature or density can be associated with 
the pixel. However, we can define the temperature and the gas density at pixels in spectra to 
be the optical depth-weighted average over the temperature and the density of all the gas that 




Fig. 1. — Temperature vs. density for random points in real space from the simulation. The 
lines are minimum absohite deviation power-law fits to each set of points in the restricted range 
1 < Ap < 2; the fitted parameters are given in Table 1. 
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contributes to absorption in the pixel. This definition will be used later to assign a gas temperature 
to any identified absorption line, which we will define to be equal to the temperature of the central 
pixel. Table 1 shows fits to the T — Ag relation for pixels, restricted over the same density range 
1 < Ag < 2: we see that the mean relation is almost the same as for random points in space. 
This is true in spite of a larger difference in the density distribution; for example, the median 
is (0.47,0.39,0.31) for random points at z = (4,3,2), and is (0.52,0.47,0.41) for spectral pixels. 
This difference in the median density is caused by the thermal broadening and velocity dispersion 
in absorbers, which spread high density regions out into low density regions. 

3. FITTING METHOD 

In this section wc describe the procedure that we use to identify and fit absorption lines. The 
temperature measurement is based on the identification of absorption lines that can be adequately 
fit by a single Gaussian in optical depth, over a certain interval around a point where the optical 
depth is maximum; the narrowest widths among the lines that can be fitted in this way will give 
us the gas temperature. Absorption systems that cannot be fitted by a single Gaussian must be 
broadened by a non-Gaussian distribution of the fluid velocity, and can therefore be discarded for 
the purpose of measuring the gas temperature. In contrast to the standard Voigt profile fitting 
approach, we make no attempt to fit the entire spectrum by superposing many absorption lines. 
Instead, we fit only small regions around minima of the transmitted flux, each one with a single 
Gaussian absorber. We therefore have a constant number of parameters to fit for each absorption 
line, making the algorithm simple, unambiguous, and fast. 

Before we explain the procedure for identifying and fitting lines, it is useful to understand 
qualitatively what the results will look like. Figure 2 shows a section of the spectrum of Q1422, 
with the transmitted flux indicated by the dotted line, and the fitting solutions indicated by the 
solid lines. Our method has selected all statistically significant maxima in optical depth that are 
well fitted by a single Gaussian and identified them as absorption lines and has discarded the rest. 
The fits are done over the regions indicated by the solid lines. We will return to this figure once 
we have described the process used to fit the lines. 

Our fitting method consists of taking each pixel in the spectrum as a candidate for containing 
the center of an absorption line. After requiring several conditions and eliminating most of the 
pixels as candidates, a final list of absorption lines is obtained, each one having a fitting window 
where a fit to the three parameters of the line (line center, central optical depth, and line width) 
is performed. None of the fitting windows from adjacent lines can overlap in the final list, as seen 
in the example of Figure 2. 

The first operation is to determine an integer window width, W, which sets the fitting region 
around each pixel P, going from P — W to P + W. The width W is the smallest one for which the 
following condition is obeyed, which ensures that there is a significant decline of the flux from the 
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Fig. 2. — The dotted line shows a piece of an observed spectrum. The sohd hnes show the fitting 
results. 



-9- 



average value at the edges of the fitting window to the hne center: 



^[F{P + W) + F{P 



W)]-F{P)>Eaa{P,P±W) , 



(2) 



where 




.1/2 



(3) 



(t{P) is the noise at pixel P, and is the first parameter of the fitting algorithm. is the number 
of "a" significance required of the flux decrease. 

There are of course some pixels where the condition in equation (2) is never obeyed for any 
width. In practice, the width W is increased only up to some maximum value Wmax before the pixel 
is discarded as a candidate for a line center. This maximum width is chosen to be large enough so 
that it does not affect the results of the algorithm. In addition, the minimum value of W is set to 
2 pixels (so that the fitting region has at least 5 pixels). Some additional parameters are used in 
the algorithm to expedite the elimination of pixels as candidates for line fits, increasing the speed 
of the code without affecting the final result; these details are described in the Appendix. 

For each pixel P where a fitting width W has been determined in this way, the line fit is 
performed by minimization. We fit the following profile to the fiux within the window: 



Here, v is the distance from the central pixel P. The three parameters of the fit are Tc (the optical 
depth at the line center), Vc (the location of the line center), and aj, (the width of the line). The 
parameter Vc is constrained to lie within the pixel P, so that when P is not close to the center of 
the line, a good fit will not be obtained. To account for the instrumental resolution, the function 
in equation (4) is convolved with a Gaussian filter of width matching the resolution of the data. 
After the fit is performed, we impose a goodness-of-fit requirement: the probability of exceeding by 
chance the value of for the best fit should be larger than a certain value Pq, which is a second 
parameter of our method. If the requirement is not satisfied, the pixel P is discarded as a candidate 
for including the center of an absorption line. The two parameters E^ and Pq can be adjusted to 
optimize the temperature measurement, based on tests using numerical simulations that we will 
present in §4. 

Because we require Vc to be within pixel P, acceptable fits to absorption lines arc only found 
in pixels that arc indeed close to a minimum of the flux, in an absorption feature that can be 
adequately fitted to a single Gaussian within a certain window. Typically, the list of pixels where 
acceptable fits are found will include groups of a few adjacent pixels around such minima. The 
final step of our algorithm is to select among any group of pixels with accepted line fits that are 
within their own fitting windows the one that yielded the best fit. This produces the final list of 
absorption lines. 
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We can now understand the fitting example shown in Figure 2. The widths of the fitting 
windows, shown by the sohd hnes, are set by the value of (here, = 12) and the noise 
level (which in this case varies from 0.004 in saturated pixels to ^ 0.01 at the continuum). 
The apparent maxima of absorption that have no corresponding fitted line do not increase in flux 
enough at their edges to satisfy the requirement in equation (2). The cluster of fitted lines near 
A = 4830 A demonstrates that our procedure does not automatically eliminate lines in blends, as 
long as they arc clearly distinct maxima. Apparently, the rcr^iiircment that the optical depth be 
consistent with a Gaussian curve is easily fulfilled by the peaks of all of the significant absorption 
lines. 

Throughout this paper, we shall be expressing all results concerning the line widths in terms 
of the equivalent temperature, denoted as B, when the line width is assumed to be due to thermal 
broadening: B = 10000 [ab / 9. 09 km s~^)^ K. This allows for an easier comparison of the results of 
numerical simulations and observations. Note that line widths have usually been presented in the 
literature in terms of the Doppler parameter, b = 2^/^(7;,. 

An important property of this algorithm is that the distribution of line widths it measures 
should converge to a fixed answer as the signal-to-noise ratio of the observations is increased. In 
the limit of negligible noise and pixel size, and perfect resolution, every true minimum of the 
transmitted flux in the spectrum should be identified, and the fitted line width should reflect 
the second derivative around the minimum, because the size of the fitting window around each 
minimum should be very small. Moreover, the second derivative around minima is a physically well 
motivated quantity to obtain the gas temperature. In contrast, the Voigt profile fitting method does 
not converge at high signal-to-noise ratio because the number of blends assumed in an absorption 
system will change. Whereas the Voigt profile method attempts to fit the entire spectrum by 
superposing lines, our new method fits only small regions around the minima of transmitted fiux 
as arising from a single absorber. 



4. APPLICATION OF THE PROFILE FITTER TO THE NUMERICAL 

SIMULATION 

4.1. Detailed Example of Fitting the Simulated Spectra 

We now apply our method to 1500 randomly selected lines of sight through the simulation 
output at z = 3, with the mean flux decrement, noise level and pixel size set to match the obser- 
vations at z = 3 (see Paper I ). We first set the two parameters of the line-fitting algorithm to 
Efi = 12 and Pq = 0.01 (we shall analyze the optimal values of these parameters in §4.2). A total of 
6378 lines are identified and successfully fitted. For each absorption line we obtain four quantities: 
the optical depth at the line center, Tc, the Doppler parameter converted to temperature units, B; 
the optical-dcpth-wcightcd temperature at the central pixel, T; and the optical-depth-weighted gas 
density at the central pixel, Ag. 



- 11 - 




0.1 1 10 100 0.1 1 10 100 



Fig. 3. — The difference, B — T, between tlie temperature estimated from the profile width and the 
true temperature at the profile center, plotted against the central optical depth Tc, for absorption 
lines in the simulated spectra at 2; = 3. Noise at the same level as the observations at z = 3 
analyzed in §6 is included in (a), and not included in (b). 
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Figure 3(a) shows B-T vs. Tc for all 6378 fitted lines. The solid lines at S - T = K and 
B — T = ±3000 K are to guide the eye in evaluating the contribution of non-thermal broadening 
to B. The noise that has been added to the simulated spectra (see §2.2.1) is responsible for most 
of the absorption lines with B < T. In Figure 3(b) we show the fitted lines from the same set of 
spectra as Figure 3(a), but without adding noise (although the noise level that each pixel should 
have was still used to weight the fits) ; the small number of lines that had 5 < T in the presence 
of noise have been almost entirely eliminated (the few remaining lines with B < T arise because the 
optical-depth-averaged temperature at the central pixel of a line can be skewed by a contribution 
from very hot gas). 

Figure 3 demonstrates that the contribution to B from fluid motions is generally large, and 
therefore the Doppler parameter of an individual line will usually overestimate the gas temperature. 
For Tc > 1, the distribution oi B — T extends to values as low as ~ 1000 K, implying that the lower 
cutoff in the B distribution should provide a good measure of the median gas temperature at a 
given Tc [assuming that the lowest gas temperatures do not extend very much below the median 
value, as is true in the simulation (Fig. 1)]. However, for weak lines (tc < 1), thermal broadening is 
never clearly dominating, preventing a model-independent estimate of the gas temperature at the 
correspondingly low gas densities. The reason why the breadth of weak lines is always dominated 
by motions is that all the low density gas has not turned around from the Hubble expansion. 

The eight histograms in Figure 4, constructed from the same fitted lines as Figure 3(a), show 
more quantitatively the cutoff on the B — T distribution near B — T = 0, for different optical 
depths. In the r ^ 0.3 and r ~ 0.6 panels of the Figure, the problem of estimating the temperature 
of the lower density gas is clearly seen. Over the range 1 < r < 20, the distribution of 5 — T is the 
desirable one for our temperature measurement: there are many lines near B — T = 0, but very 
few with B <T. 

The sharpness of the cutoff in Figure 4 should of course be degraded in the observable B 
distribution, owing to the scatter in the temperatures of the lines. Figure 5 shows histograms of 
B along with histograms of T. By observing the histogram of B we would like to determine the 
median of T. For Tc > 1, the cutoff on the distribution of B appears to coincide well with the peak 
of the distribution of T. As we saw clearly in Figure 3, the number of lines with B c^T decreases 
quickly for Tc < 1. 

It is interesting to note in Figure 5 the change in the distribution of the true temperature 
in the simulation with Tc- At low optical depth, the T distribution is very narrow as a result of 
the simplicity of the evolution history of gas at low density, which generally expands peacefully 
in the voids, heated by photoionization and cooled adiabatically by its expansion. As the density 
increases, the heating history of the gas becomes more heterogeneous. The gas at higher densities 
is shock-heated more frequently and to a greater degree, and the evolution of the density itself since 
reionization (when the initial temperature is set) is more highly variable. For r > 10 the range 
of gas temperatures increases, making the interpretation of the cutoff on the B distribution more 
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Fig. 4. — Distribution of the non-thermal broadening contribution to the Hne widths, B — T, at 
eight bins of the optical depth at the line center, Tc- N is the total number of lines in each Tc bin. 
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Fig. 5. — The solid lines show the histograms of B and the dotted lines show the histograms of T. 
The fitted lines used here are the same ones shown in the previous several Figures (described in 
§4.1). The probabilities in the T histograms are all reduced by a factor of 2.5 to facilitate visual 
comparisons. 
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ambiguous. 

4.2. Optimizing the Fitting Control Parameters 

The two parameters of the hne fitter, and Pq, should be optimized to give the best statistical 
error bars on the location of the lower cutoff of the distribution of B, which we will use to estimate 
the gas temperature. In order to do this optimization, we define a quality measure, Q = (Ng — 
Nb){Ng + Nh)-^/^, where Ng and Nb are the number of fitted hues with OK < B-T < 3000 K, and 
with —3000 K<5 — r<OK, respectively. Basically, Q is a measure of the statistical significance 
of the increase in the number of lines as the B = T cutoff is crossed, computed by comparing the 
number of lines in bins of width 3000K on each side of the cutoff. The larger the value of Q, the 
more accurately we should be able to determine the temperature. For example, the set of fitted 
lines shown in Figures 3-5 (fitted using = 12 and Pq = 0.01) gives Ng = 281 and Ni, = 40, 
resulting in Q = 13.5. When we remove the noise in Figure 3(b) we find Ng = 236 and N^ = 1, 
yielding Q = 15.3. We use a bin width of 3000 K because the errors in the measured temperature 
from our data will be of about this magnitude (see §6). 

Table 2 shows the Q values for a broad range of values of E^, and Pq = 0.01. We also list 
the Q value using two different values of Pq, removing the noise from the spectra, removing the 
continuum fitting approximation, and using different random seeds for the added noise. We find 
that the changes in Q are usually smaller than the changes that can result from simply using a 
different set of random numbers for the noise that is added to the spectra. We conclude that the 
precise values taken by the parameters are not actually very important to the results. We use 
Ed = 12 and Pq = 0.01 when we analyze spectra matching the z = 3 observational properties in 
the rest of this paper. Prom the range of E^ with Q close to its maximum, we chose the smallest 
value of Ed, Ed = 12, because we expect that increasing the number of accepted lines will make 
the analysis procedure more robust, particularly the computation of the error bars. 

We use two other redshift bins for the observational data (see §2.1), z = 2.4 and z = 3.9. 

Because the mean flux decrement, pixel size, and noise level are different in each bin, we determine 
a best value of Ed separately for each. Tests similar to the one in Table 2, using the mean flux 
decrement, noise level, and pixel size matching the observations in the high and low z bins, show 
that Ed = 9 and Ed = 8 are the best values to use when analyzing data in the z = 2.4 and z = 3.9 
bins, respectively, although the values of Q obtained are only weakly sensitive to Ed in each case. 
We fix Pq = 0.01 at all three redshifts, because changing it does not significantly increase Q. 
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5. TESTS OF THE DERIVATION OF THE IGM TEMPERATURE USING THE 

SIMULATION 

In this section we combine our line fitting method with elements of the method of Bryan &; 
Machacek (2000) for constraining the temperature of the IGM by measuring the lower cutoff on the 
distribution of B. First we define the cutoff and how we associate it with an estimated temperature. 
Then we use the simulation to translate the observed temperature-optical depth relation into the 
desired T — Ag relation . All the results in this Section are obtained from the simulation. In §6 
we present results for the temperature measured from the observational data, using the method 
described in this Section. 



5.1. Locating the Cutoff on the Distribution of B 

Bryan &; Machacek (2000) presented a useful technique for quantifying the lower cutoff on a 5 
histogram like those in Figure 5. They smooth the histogram with a Gaussian filter and define the 
cutoff to be the location of the maximum of the first derivative of the smoothed histogram. This 
derivative is given by 



dP ^ 

— ^ (x^(Sj -Bi)exp 



l{B^jf 
2 al 



(5) 



where the fth bin has temperature Sj, j is the label for N individual fitted absorption lines with 
fitted temperature and as is the smoothing length. We plot in Figure 6 the smooth histogram 
of B and its derivative, as well as the smoothed histogram of the gas temperature, for gb = 5000 
K, using spectra from the z = 3 simulation output with the z = 3 observational noise, pixel size 
and mean flux decrement. 

We define the Doppler parameter cutoff, Bc{tc), to be the value of B where dP/dB, given 
by equation (5), is maximum, for the optical depth bin labeled by Tc- Our estimate of the gas 
temperature at optical depth Tc is Bc{tc), after applying a small correction that we describe in 
detail in the remainder of this section. 



5.1.1. Error Bars on Bc{tc) 

We compute error bars on the location of the cutoff on the B distribution {Be) by bootstrap 
resampling (Press ct al. 1992). We generate a bootstrap realization of the B histogram, for an 
optical depth bin containing fitted lines, by selecting N lines at random from those in the bin 
(with replacement) and recomputing the histogram from the new set of lines. The error on Bq 
is given by the dispersion in the Be values measured from many bootstrap realizations of the 
histogram. 
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Fig. 6. — The derivative of the smoothed histogram of B (solid line), the smoothed histogram of 
B {dashed line), and the smoothed histogram of T {dotted line), from the absorption hnes fitted 
in the z = 3 simulation output with the fitting control parameters = 12 and Pq = 0.01. The 
smoothing length is as = 5000 K. The curves are arbitrarily normalized for clarity. 
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5.2. Comparison Between Bc{tc) and the Temperatures in the Simulation 

In this subsection we investigate the relationship between Bc{tc) and the physical temperature, 
T, of the absorbers with central optical depth Tc- We use spectra from the simulation, where 
we know the optical-depth-weighted temperature and density, T and A, of the absorption lines. 
Wc match the simulated spectra to the mean flux decrement, noise level and pixel size of the 
observations at z = 3 (see §2.2.1). In §5.3 we show how the comparisons change when the simulated 
spectra are matched to the properties of the z = 3.9 and z = 2.4 observations. 

Figure 7 explores the meaning of Bc{tc) in detail, using the z = 3 simulation results. We 
have separated all the fitted lines into 10 bins of the central optical depth, Tc, choosing the bins 
so that each contains an equal number of absorption lines. Wc smooth the B histogram with 
as = 5000 km s~^ to compute the cutoff temperature in every optical depth bin. Then, we compute 
the median central optical depth and gas temperature of the set of absorption lines that satisfy 
\B — Bc\ < 5000K, in each optical depth bin. These sets of lines are the ones that actually determine 
the location of the cutoff in the B histogram, so we examine first the relation between Be and 
their median temperature, which we denote by Tp. The values of Bq are shown as crosses with 
error bars. ^ The temperature Tp is shown as the filled squares (the errors on the temperature are 
much smaller than those on Be)- 

The triangles in Figure 7 show the effect on the derived cutoff. Be, of reducing the histogram 
smoothing to ap = 3000 K. There is no important difference with the crosses (the error bars for 
this smaller smoothing are similar to those on the crosses). By experimenting with different values 
of ap, we have found that Be is not affected if ap is reduced, although the error bars obtained are 

significantly larger when reducing it below crp = 3000 km s^^. Increasing ap beyond 5000 km s~^ 
leads to an increase of Bq, because the smoothing is then larger than the intrinsic sharpness of the 
cutoff in the B distribution. We therefore adopt crp = 5000K from this point forward in the paper. 

The open squares show the effect of using all the absorption lines in every optical depth bin to 
compute the median temperature and optical depth, rather than using lines with \B — Bc\ < 5000K 
only. There is a negligible difference in the derived temperature of the absorbers as a function of Tc, 
the reason is that, as we have previously seen, broader lines are by and large systems with higher 
fluid velocity dispersions, but their gas temperatures are not significantly greater, except at the 
highest optical depths where there is a slight difference (the systematic shift to the right of the 
open squares relative to the filled ones is due to a larger median optical depth of the broad lines 
within each optical depth bin). In the rest of the paper, we always compute the medians of any 
properties of the absorption lines using only lines with \B — Bc\ < 5000K. 



*The error bars are from bootstrap realizations on 1500 lines of sight. In reality, these error bars may be slightly 
underestimated because the mean separation between 1500 hncs of sight in the simulation wc use is only ~ 60kms~^, 
comparable to the flux correlation length in the spectra (Paper I ). Obtaining better statistics of the theoretical 
prediction for Be would require a larger simulation. However, for the analysis in the present paper the errors in the 
determination of the temperature from the observations are much larger. 
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Fig. 7. — The crosses with error bars show Bc{tc), computed from the simulated spectra at z = 3, 
and using as = 5000 km s~^. The open triangles show Be when the smoothing of the B histogram 
is reduced to ctb = 3000 K. The open squares are the median T of all the absorption lines in each 
bin, while the filled squares are the median T of lines with \B — Bc\ < 5000K only. The dashed 
line shows Tr, the median value of T for random cells in the simulation, and the dotted lines show 
the 25th and 75th percentiles. 
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So far, we have seen that the Doppler parameter cutoff Be provides a good estimator for the 
gas temperature of absorption systems at a given optical depth. Our ultimate goal, however, is to 
measure the median gas temperature at a given gas density, for randomly selected points in the 
IGM, which we shall henceforth refer to as Tr. The peaks in absorption are at special locations, 
so their median temperature, Tp, will generally not be exactly the same as Tr. To examine this 
question, we first compute the median density of the absorption lines in each optical depth bin, 
and then we calculate the median temperature of randomly selected points at this gas density. The 
result is shown as the dashed line in Figure 7; the dotted lines give the 25 and 75 percentiles of 
the temperature distribution at random points with the same gas density. Comparing the dashed 
line to the filled squares, we see that the fitted absorption lines with relatively large optical depth 
typically have Tp < Tr, i.e., they are colder than random fluid elements at the same density. We 
believe the reason for this effect is the characteristic double-shock structure around the absorbers 
(Cen et al. 1994): the gas in the highest density tube along a filament (or the highest density 
surface along a sheet) is located between shocks, so it has been subject to less shock-heating than 
the surrounding gas. 

Our method of analysis of the observational data in §6 will automatically correct for this differ- 
ence between the temperature of the absorption lines and the temperature at random points. This 
systematic difference may introduce a potential uncertainty in the derivation of the gas tempera- 
ture if it depends on quantities like the resolution of the simulation, the cosmological model that 
is assumed, or the heating at the reionization epoch. However, the temperature difference between 
lines and random points is negligible compared to the observational errors we will compute for the 
temperature in §6, at least in the simulation we analyze here. 

5.2.1. Testing With Other T-Ag Relations 

We need to test the robustness of the finding that Bc{tc) — Tp{tc) for Tc > 1, which we have 
established so far in tests on the z = 3 simulation output. We can change the T — Ag relation that 
we are measuring by simply using the z = 2 oi z = A outputs from the simulation (see Table 1), 
still creating spectra with mean flux, noise level, and pixel size matching the z = 3 observations, 
as described in §2.2.1. In addition to the varying T — Ag relations, these spectra have different 
amplitudes of fluctuations, and different Hubble constants. 

Figure 8 shows the comparison between Bc{tc) and Tp{tc), using the z = 2 and z = A outputs 
of the simulation in addition to z = 3. The {pentagons, squares, triangles) show Be for the z =(4, 
3, 2) simulation output, and the (solid, long-dashed, short-dashed) line show Tp{tc)- We see that 
Be traces the temperature changes at the different redshifts extremely well, tracking the different 
slopes aX z = 2 and z = 4 perfectly, and matching the increased overall temperature at z = 3. 
All three redshift outputs show the same strong increase in Be above the actual temperature for 
optical depths below unity, corresponding to Ag ~ 1 in each case. 
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Fig. 8. — (Pentagons, squares, triangles): Be for the z =(4, 3, 2) simulation output. (Solid, 
long-dashed, short-dashed) lines: Tp (median temperature of fitted lines) for z =(4, 3, 2). All 
simulated spectra have the mean flux decrement, pixel size and noise level set to match the z = 3.0 
observations. 
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5.3. The Bc-T Relation at Different Redshifts 

We have seen that Be traces Tp remarkably well for the mean flux decrement and noise level of 
our observational data at z = 3. We shall now verify that this is also true when the flux decrement 
and noise is set to the values appropriate for the other two redshift bins in which we separate the 
data, at z = 2.41 and z = 3.89. 

We first introduce a new type of figure that shows more clearly how accurately Bc{tc) traces 
the gas temperature. Figure 9(a) shows Be — Tr, where Tr is the median temperature of random 
points at the density of each optical depth bin (recall that this density is defined as the median 
density of absorption lines in each optical depth bin that have \B — Bc\ < 5000K). The (pentagons, 
squares, triangles) are obtained from the simulation outputs at z =(4, 3, 2), and are the same 
points as in Figure 8 except that we have subtracted Tr. The temperature is correctly traced by 
Be at Tc > 1, and drops below Be at lower optical depths in the same way, independently of which 
simulation output we use. At Tc > 10, Be falls significantly below the temperature in the z = 2 
output (this is seen clearly using more optical depth bins). The reason is that the temperature 
dispersion is higher at z = 2 in the simulation, causing Be to reflect the lower cutoff of the true 
temperature distribution. 

The results when we fix the mean flux decrement and noise to the z = 2.41 redshift bin of the 
observational data (taken from Paper I ) are shown in Figure 9(b). Here, we use the parameters 
Ed = 9 and Pq = 0.01 for the line fitting algorithm (see §3). The results are shown also using 
all three simulation outputs. Again, we find the temperature is well traced by Be, but over a 
range of Tc that has shifted to lower values. This is a result of the decreased optical depth at 
fixed gas density. Absorption systems that are primarily thermally broadened exist above a gas 
density that is approximately constant at each redshift, but the corresponding optical depth varies 
rapidly with redshift. At high optical depths, an additional effect is important in changing the 
degree to which Be traces the gas temperature: the increased shock-heating at low redshifts (with 
increasing velocities of collapse) implies a higher temperature dispersion, even at a fixed gas density. 
Therefore, Be drops further below the median gas temperature as the redshift decreases. 

The result for z = 3.9 is shown in Figure 9(c) (we use = 8 at this redshift). The Doppler 
parameter cutoff (Be) now traces the temperature only at high optical depth, Tc > 5, for exactly 
the same reason: a fixed gas density has shifted to a significantly higher optical depth due to the 
increase in the mean flux decrement. 

5.4. Determining Tji{Ag) Using Be{Tc) 

So far, we have seen that the Doppler parameter cutoff traces the gas temperature over a 
reasonable range of optical depth. We have shown this to be a consequence of the presence of some 
absorption lines that are primarily thermally broadened, and of the small dispersion of the gas 
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Fig. 9. — {Pentagons, squares, triangles): Be — Tr, at the z =(4, 3, 2) simulation output. The 
mean flux decrement, pixel size, and noise level is set to match the z =(3.0, 2.4, 3.9) observations 
in panels (a, b, c). 
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temperature. We therefore expect that this relation between Be and Tr will not be significantly 
changed depending on the model adopted in the simulation, or when the numerical resolution is 
increased. 

However, even over a restricted range of optical depth where Be and Tr are best matched, 
the simulation predicts a difference between them, which we want to correct for when analyzing 
the observational data (although this correction could be model dependent, and will need to be 
compared with other simulations in future work). In addition, we need to relate the central optical 
depth Tc to the optical-depth- weighted gas density of the absorber. A;,, in order to derive the T — Ag 
relation of the gas from the observed Be as a function of Tc- Given the limited amount of data 
that we will analyze in this paper, it will be sufficient to parameterize the T — Ag relation by a 
power-law, 

T = mAg/A,y-' . (6) 

where is chosen to make the error on T^, and 7 — 1 uncorrelatcd. As discussed in §2.2.2, this 
power-law should be understood only as a fit to the true relation, which should be more complex. 
The power-law form will be adequate here given our error bars, but larger sets of data might be 
used to detect deviations from a power-law. 

In this subsection we develop the method to derive the parameters and 7, given the deter- 
mination of Be at different optical depths. 

5.4- 1- Accounting for the Systematic Offsets Be — Tr 

To obtain an accurate estimate of the gas temperature from observations, the systematic 
differences between Be and Tr shown in Figure 9 can be used to correct the observed Be- However, 
this correction may be dependent on the model, and this dependence will not be known until a 
wide variety of additional simulations are analyzed. We therefore use only absorption lines over 
the range of Tc where the offset between Be and Tr is small {\Be — Tji\ < 5000K for all three 
simulation outputs). The following optical depth ranges will be used: 0.41 < Tc < 5.4 for z = 2.4, 
1.0 < Tc < 19 for z = 3.0, and 3.8 < Tc < 47 for z = 3.9. 

For each fitted line obtained from the observations (within the accepted optical depth range), 
we determine a temperature correction at the optical depth of the line by linearly interpolating 
from the two adjacent points in Figure 9. The set of points used depends on the redshift bin of 
the observations and the simulation output we choose for the analysis. The corrected line width is 
B' = B — AT{tc), where B is the observed line width, AT{tc) = Bc,s{'Tc) ^ T^^siTc), Bc.si'Tc) is 
the cutoff of the B distribution in the simulation, and Tji^s{tc) is the median temperature in the 
simulation at random points with gas density equal to the median density of the absorption lines 
satisfying \B - Be{Tc)\ < 5000K. 
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5.4-2. Translating Tc into Ag 

After the corrections just described, we have an estimate of the temperature as a function 
of Tc. What we need, however, is the temperature as a function of A^. This will be derived by 
using a transformation from Tc to Ag that we obtain from the simulation. This introduces an 
inevitable model dependence in our mcasTirement: the relation between Tc and A^, given a fixed 
flux decrement, should essentially be subject to the same uncertainties that appear in deriving the 
parameter fi ^ (Qf^h'^)'^ /T / H(z) (where T is the photoionization rate due to the cosmic background) 
from the observed mean flux decrement (see Ranch et al. 1997, Paper I ). 

Figure 10 is a scatter plot of the density and optical depth of the lines fitted from the z = 3 
simulation output (with mean flux decrement matching the z = 3 observations). The crosses show 
Ag vs. Tc for all of the lines, while the filled squares show only lines satisfying \B — Bc{tc)\ < 5000K, 
where Bc{tc) was determined for 10 optical depth bins as described earlier. The lines that determine 
the B cutoff (shown as squares) tend to have a higher optical depth than other lines at the same 
density, because of their lower velocity dispersion. 

To obtain the Ag — Tc relation, we assign the median density of the lines in each optical depth 
bin that wc use to measure the temperature (i.e., those with \B — Bc{tc)\ < 5000K) with the 
median Tc of the same lines. We use interpolation to calculate the density corresponding to any 
value of Tc for the fitted lines in the observed spectra. 

One of the quantities affecting the Tc — A^ relationship is the density-temperature relation 
itself, essentially because the temperature affects the recombination coefficient, which then changes 
the neutral fraction at a given density. In order for our determination of the T — Ag relation to be 
self-consistent, we need to change the temperatures in the simulation so that they agree with the 
same T — Ag relation . We do this in the following way: after determining a preliminary T — Ag 
relation using the simulation with the true temperatures, we modify the temperature in every cell 
of the simulation using the formula 

Tl = Ti- UAi/A,y-' + r^(A,/A,)(^-')' , (7) 

where Tj is the original temperature at cell i, T/ is the modified temperature, and 7 — 1 are the 
parameters of the original T — Ag relation of the simulation, and and (7 — 1)' are the parameters 
of the new T — Ag relation that we wish to impose. We then iterate the application of this formula 
until the modified T — Ag relation of the simulation matches the one from the observations. This 
modification of the temperatures in the simulation causes only a small change in the derived T — Ag 
relation (the value of Tq is modified by only ~ 5%) . 

5.4.3. Fitting for 7 — 1 Without Binning 

We now describe the method we use to fit the parameters T^, and 7 — 1 to the values of A^ 
and B' of a set of fitted lines. The simplest method would be to separate the lines in density bins, 




Fig. 10. — Ag vs. Tc for a set of fitted absorption lines from the z = 3 simulation output. Crosses 
mark all the lines, squares indicate the lines with \B — Bc{tc)\ < 5000K. 
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measure the cutoff Be in every bin, and tfien fit the power-law relation to the values of Be obtained 
at every bin. However, the binning could result in a degradation of the measurement errors: at 
least 50 lines arc needed to obtain a reasonable estimate of the cutoff Be, and since our data yield 
only a few hundred lines for each redshift bin, the binning in line density would need to be very 
coarse. 

There is a simple solution to this binning problem if the cutoff on the distribution of fitted 
lines, in the B' — Ag plane, is described by a power-law (recall that B' is the width of each fitted 
profile minus the expected non-thermal broadening at its optical depth, as described in §5.4.1). 
After we have associated a gas density with each fitted line and corrected their temperatures 
using the simulation predictions, we rotate the absorption lines in the B' — Ag plane for many 
different assumed values of 7 — 1, using the formula 

B" = B'- n{Ag/A,r-' . (8) 

For each assumed 7— 1, we apply the cutoff determination technique to the B" distribution (without 
density binning) to find a value for the temperature and a maximum value for dP/dB [see eq. (5)]. 
As 7 — 1 is varied, the best fit value of 7 — 1 is the one that results in the maximum value of dP/dB, 
i.e., the sharpest cutoff on the B" distribution. 

Note that the value of T^, used in equation (8) affects the size of the temperature changes in the 
rotation. We therefore also iterate in the determination of and 7 by this procedure. In practice, 
the measurement of is barely affected by the rotation in equation (8), so a single iteration is 
sufficient. 

Before presenting the results of applying our method to the observations, it will be useful at 
this point to summarize the full procedure we have described for measuring and 7 — 1 from the 
B distribution of the fitted absorption lines at every redshift bin. This consists of the following 
steps: 

1. Eliminate the absorption lines with Tc outside the range where the temperature measurement is 
expected to be effective. 

2. Correct the values of B for all of the remaining lines using the systematic offset predicted by 

the simulation (Figure 9). 

3. Use the Tc — Ag relation in the simulation to associate a value of Ag with each fitted line. 

4. Determine an initial estimate for T^, and 7 — 1 with the method of fitting the power-law cutoff 
in the B distribution as a function of Tc that was just described. 

5. Modify the T — Ag relation in the simulation to more closely approximate the relation measured 
from the observations. 

6. Repeat steps 2-5 until the T — Ag relation measured from the observations matches the relation 
in the simulation. 
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6. ANALYSIS OF THE OBSERVED SPECTRA 

We now apply our line-fitting method to the observed spectra. Table 3 lists, for each redshift 
bin, the redshift range (^^^ri^n to Z'jyiax) and mean redshift the mean flux decrement the mean 
noise, the total number of pixels, the total path length, the value of we use to fit the lines, 
the total number of lines fitted, and the number of lines that we actually use to measure the gas 
temperature in the range of central optical depth from Tmin to Tmax- 

The parameters {B,Tc) of all the fitted lines axe shown in Figure ll(a,b,c) for the redshift bins 
z =(3.9, 3.0, 2.4), as outlined crosses when the lines are not in any of the regions suspected to 
include metal lines, and as simple crosses when they are. For the quasar KP 77 (included in the 
redshift bin z = 2.4), the analysis to identify potential metal lines was not done, so all lines from 
this quasar are shown as outlined crosses. The lower cutoff on the B distribution is clearly visible 
to the eye at all three redshifts, especially when the metal lines are ignored. 

In Figure 12(a,b,c) we compare the B histograms of the observed lines, within the ranges of 
optical depth that we will use for the temperature measurements, to the B histograms of fitted 
lines from the simulation outputs with redshifts closest to the means of the observations, in the 
same optical depth ranges. The observed absorption lines obviously have higher temperatures than 
the simulated ones in all three cases. 

It is interesting to compare directly Bc{tc) from the observations and from the simulation, 
before wc determine T = T^(A(,/A^)'''~^ using the more involved method described in §5.4. Figure 
13(a,b,c) shows the values and errors of Be measured from the observations and the simulation 
output that is nearest in redshift, using lines over the optical depth bins indicated by the horizontal 
error bars (the optical depth bins contain equal numbers of observed lines). The vertical dotted 
lines indicate the optical depth range that we use for the final temperature measurement (for the 
z = 3.9 analysis, the upper limit on is outside of the figure, and eliminates a negligible number of 
lines). Recall (Figure 9) that we do not expect the points with lower Tc to accurately reflect the real 
temperature, except at the lowest redshift. The observed lines again appear to be hotter than the 
simulated lines. These results are listed in Table 4, along with the temperature offset AT used as 
a correction to the temperature (see §5.4.1) and the estimated gas densities, Ag, that we find once 
the T — Ag relation in the simulation has been adjusted to match the observed one (determined 
below). If more observed spectra were available, this method of binning in optical depth would be 
preferable to the method in §5.4 because it does not require the assumption of a power-law T — Ag 
relation . Each bin in optical depth would be associated with the density listed in Table 4, and the 
value of Be would be corrected by the listed AT. 

We now determine and 7 — 1 by the method described in §5.4. In order to obtain the 
T — Ag relation and the temperature offsets AT, we can use any of the three simulation outputs 
at 2; = 2,3,4 for any of the three redshift bins in which we have divided the data (although we 
change the mean flux decrement of the simulated spectra to the observed one at each redshift 
bin). The different redshift outputs of the simulation are approximately equivalent to assuming 




Fig. 11. — (a, b, c) Fitted lines in the observed spectra at z =(3.9, 3.0, 2.4). Outlined crosses are 
lines that are not in regions containing potential metal lines. 
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Fig. 12. — Distribution of B of fitted absorption lines, (a) Solid line: observed lines at z = 3.9, 
with 3.8 < Tc < 47. Dotted line: Lines from simulation output at z = 4, with 3.8 < Tc < 47. (b) 

Same as (a) for the z = 3 redshift bin and z = 3 simulation output, with 1.0 < Tc < 19. (c) Same 
as (a) for the z = 2.4 redshift bin and z = 2 simulation output, with 0.41 < Tc < 5.4. Nq is the 
total number of observed lines in each histogram. 
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Fig. 13. — (a, b, c) Crosses with large error bars: Bc{tc) for observed lines at z =(3.9, 3.0, 2.4). 
Squares with small error bars: Bc{tc) from the simulation at z =(4, 3, 2). Only lines in the range 
of Tc between the vertical dotted lines [or to the right of the single dotted line in (a)] are used for 
our final temperature measurement. 
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different models with a different amplitude of the power spectrum (Paper I ), so we can use the 

different outputs to check that our measurement of the temperature is not greatly sensitive to the 
model that is assumed. Paper I found that the amplitude of the initial density perturbations in our 
simulation needs to be reduced by about 15% to agree with the observed power spectrum of the 
transmitted flux, meaning that the simulation output at z = 4 has fluctuations that most closely 
match the observed ones at z = 3, and are slightly higher than the observed fluctuations at z = 3.9. 
The z = 2.4 observational bin is closest in amplitude to the simulation output at z = 3. The 
most reliable temperature results should therefore be obtained by using the z =(4, 4, 3) simulation 
outputs to analyze the z =(3.9, 3.0, 2.4) observations, but we shall also give results for z =(3.0, 
2.4) analyzed using the z =(3, 2) simulation outputs. 

The result at z = 3, using the z = 4 simulation output to predict the r — Ag relation and the 
necessary correction for non-thermal broadening, is: 

/ A \ [0.29±0.30] 

T = [20200 ± 1300] K , (9) 

where the error bars on and 7 — 1 arc uncorrelated (which defines A^). The error bar on 
reflects only the uncertainty in the Paper I determination of the mean flux decrement, which 
affects the relation between density and optical depth. This result implies Tq = 18400 ± 2100K, 
where Tq is the temperature at the mean density. The pivot density, A,^ = 1-37, corresponds to 
an optical depth = 1.83 (from the relation obtained as described in §5.4.2) When we repeat the 
fltting using the z = 3 simulation output we flnd T = [19600±1500](Ag/[1.24±0.10])[°-'^-'^=^°-^^lK, or 
Tq = 18300ib 1800K (t+ = 1.74). The difference between the two values for A^, 1.37 using the z = 4 
simulation, and 1.24 using the z = 3 simulation, are mostly a reflection of the different optical depth 
normalizing factors (i.e., rescalings of the baryon density or the strength of the ionizing background) 
needed to match the observed mean flux decrement. The normalizing factor is smaller for the z = 4 
simulation, giving a larger for the same optical depth, because the density fluctuations are of 
lower amplitude, leading to less saturated absorption and more absorption in voids (see Paper I for 
a more detailed discussion of the optical depth normalizing factor). 

All the results obtained at the three rcdshift bins are listed in Table 5. The data analysis at 
z = 2.4 and z = 3.9 is similar to the analysis at z = 3, except for the differences that we mention 
below. 

The temperature results at z = 2.4 differ by ~ 2000 K when the z = 2 simulation output is 
used instead of the z = 3 output, once the measured values of for the two are extrapolated to 
the same density. Most of this difference results from the difference in A^, and the relatively high 
value of 7 — 1 that is obtained at z = 2.4. The difference in A,^ is caused by the different amplitudes 
of density fluctuations in the two simulation outputs, so we expect that the temperature derived 
using the z = 3 simulation output, which has the correct amplitude of fluctuations, is more reliable 
(see Paper I ). 

In order to see the evolution of the temperature with redshift, we need to obtain the temper- 
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ature at a fixed over-density at each redshift. It is useful to obtain the temperature at the mean 
density, Tq, to compare our resuhs to other work. However, the vahics of are close to = 1.4 
at all three redshift bins, and we can therefore have a more robust result for the temperature 
evolution if we examine the temperature at A^ = 1.4, which we denote as T1.4 in Table 5. 

We note here that, because of the small range of optical depth used at z = 3.9, at this redshift 
we were forced to smooth the B" histograms with a Gaussian filter of width aB = 7000K, instead 
of our standard as = 5000K, in order to avoid problems with multiple, approximately equivalent, 
maxima of equation (5) as 7 — 1 is varied. 

The primary results of this paper, the measurements of ri.4 and 7 — 1, are summarized in 
Figures 14 and 15, respectively. These show two important conclusions: first, the temperatures are 
higher than the value expected if photoionization heating in equilibrium is the only heating source. 
Second, we find no evidence for a rapid change of the temperature with redshift. 

7. DISCUSSION 

This paper presents a measurement of the temperature-density relation of the intergalactic gas 
in the redshift range 2.4 < z < 4. The new method we have developed to perform this measurement 
is based on the same general idea as the previous work by Schaye et al. (1999), Ricotti et al. (2000), 
and Bryan Sz Machacek (2000): provided that there is a tight relation between the temperature 
and density of the gas, absorption lines of similar central optical depth should have little dispersion 
in their thermal broadening, and the varying line widths should correspond to variable amounts 
of hydrodynamic broadening. Occasionally, some absorption lines will be subject to only a small 
degree of hydrodynamic broadening; this will typically happen when most of the atomic hydrogen 
occurs near a velocity caustic along the line of sight. We therefore expect the histogram of line 
widths to show a rapid increase near the value of the Doppler parameter corresponding to the 
gas temperature. In the absence of noise, every line should be wider than the thermal broadening 
width, so at least an upper limit to the temperature can be obtained unambiguously. 

The tests we have performed using a numerical simulation of the Lya forest, based on a CDM 
model that successfully reproduces the observations of large-scale structure at present, confirm this 
general idea. However, they show that this method to recover the gas temperature works efficiently 
only over a limited range of line optical depths, which corresponds approximately to a range of gas 
over-density 1 < A^ < 3. At lower densities, the gas is generally in Hubble expansion and this 
effect dominates the contribution to the line widths in essentially all the lines. The minimum line 
widths can therefore only provide an upper limit to the temperature of this low-density gas. Of 
course, the simulation can in principle be used to correct for the effect of line broadening due to 
Hubble expansion, and to obtain the temperature by subtracting the hydrodynamic contribution 
to the minimum line widths. However, the results can then strongly depend on the assumed model 
and the numerical resolution of the simulation, especially as the thermal broadening becomes a 



2.5x10* 



2x10" - 



1.5x10* - 



10* - 



5000 - 



Q I I I I I I I I I \ I \ I I I I I I I I I I I I I 

4.5 4 3.5 3 2.5 2 



Fig. 14. — The observed temperature at = 1.4. The two points at z = 2.4 and z = 3.0 (offset 
slightly to distinguish them) show the result of using two different simulation outputs to analyze 
each of these redshift bins. The left and right points at z = 2.4 are for the z = 3 and z = 2 outputs, 
respectively. The left and right points at z = 3.0 arc for the z = 4 and z = 3 simulation outputs. 
The vertical dotted lines show the boundaries of the redshift bins. 
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Similar to Figure 14 except here we plot the results for 7 — 1. 
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small effect compared to the expansion. 

At very high densities, the increasing dispersion of the temperature at a given gas density can 
result in a large difference between the typical gas temperature and the "lower cutoff" in the line 
width histogram. This implies again that the recovery of the median gas temperature from the 
distribution of line widths is highly sensitive to the model assumptions that affect the temperature 
dispersion and the turbulent motions in the gas. 

There are two main differences between the method we use to measure the gas temperature, 
and that used by previous authors. First, our line detection algorithm avoids the necessity of the 
Voigt-profile fitting method to "deblend" lines, by simply throwing out any "absorption lines" that 
do not correspond to a clearly identified minimum in the transmitted flux, or that do not have a 
large enough region around that minimum that is adequately fitted by a simple Gaussian in optical 
depth. In the method we use here, every line width is essentially a measurement of the second 
derivative around a minimum in the flux. This is one important reason why the total number of 
lines we identify is much lower than Schaye et al. (2000), even though we use nearly the same data 
set. The second difference is that we restrict the absorption lines we use to lie within the range of 
central optical depth where the correction that needs to be applied to the temperature measured 
from the histogram of line widths (as described in §5.4.1) is small. 

These differences explain our substantially increased error bars in measuring Tq and 7, relative 
to those of Schaye et al. (2000). However, we believe our error bars are more reliable and model- 
independent, for the reasons we have discussed. We also note here that, even though we cannot 
rule out a substantial change of the temperature we have measured depending on the numerical 
simulation of the Lya forest that is used for comparing to the observational data, this possibility 
appears unlikely for several reasons, in addition to the arguments explained before about the 
small size of the correction AT that needs to be applied to the line width cutoff Be to obtain 
the gas temperature (see §5.4.1). The gas temperature of the simulation we use is lower than 
that observed only by a small amount: in the simulation, the temperature at the mean density is 
ro=(14000, 15800, 12800)K at z =(4, 3, 2), while our measurement from the observational data is 
To = (17400 ±3900, 18400 ±2100, 17400 ± 1900)K for z =(3.9, 3.0, 2.4). The amplitude of the flux 
power spectrum of the simulation is also very close to the observationally determined one (Paper 
I ). We have also shown that we obtain nearly identical results for the temperature measurement 
when using different simulation outputs to analyze the same observations (see Figures 14 and 15) 
meaning that the dependence on the amplitude of the power spectrum is weak. Our result for the 
gas temperature might also be affected by the limited resolution of the simulation we use (with 
a comoving cell size of 35 Kpc). We have not yet performed a convergence test for the effects of 
resolution on the Lya forest; however, Schaye et al. (2000) find that a mean particle spacing of 
~ 45 Kpc in their SPH simulations is sufficient for convergence. 

We now compare our results for the evolution of the gas temperature with previous measure- 
ments. Ricotti et al. (2000) find a temperature at the mean density (from their Figure 12b) 
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To ~( 186001^^^^°, 23400^^^^^°, 17000l^^^g°)K, at z ~ (3.6, 2.75, 1.9). Considering their large 
statistical error bars, our results appear to agree well with theirs, although the true temperature 
at z ~ 2.75 must be at the low end of their error bar. Schaye et al. (2000) give their temperature 
results as 16 separate points, in their Figure 5. We read off their points in each of our redshift 
bins and create error weighted averages for comparison with our results, finding Tq ~(12200ib 1700, 
17300±1400, 14000±1300) at z ~(3.75, 3.2, 2.4). In order to compare these with our temperatures, 
which have lower errors at = 1.4, we extrapolate the Schaye et al. (2000) results to = 1.4 
using their measured value of 7 — 1, obtaining T1.4 ~ 13600 it 2000K z ^ 3.75. Our result at 
z = 3.9, Ti,4 = 20100 ± 2800K, is higher than theirs by ~ 1.9a. Because of this, we do not find 
evidence for the increase of the temperature with time between z = A and z = 3 that Schaye et al. 
(2000) reported. Our results are consistent with a constant temperature. 

Actually, the data set analyzed by Schaye et al. (2000) is almost identical to the one we 
analyze in this paper, with 7 of the 10 quasars used in the two papers being identical. However, our 
methods of analysis are very different, so that even the statistical error bars from the two analysis 
are largely independent. The main difference, as mentioned before, is in the number of lines that 
are identified. For example, we use 98 lines to measure the temperature in our z > 3.4 bin, while 
Schaye et al. (2000) use about 550 lines in a comparable bin, even though ~ 80% of the data that 
they use is identical to ours. 



7.1. Hell Reionization as a Heating Mechanism 

Is the value of the temperature we have measured in agreement with the known sources of 
heating and cooling in the intergalactic gas? The evolution of the temperature is determined by 
the equation: 

where the cooling and heating rates per particle are denoted as follows: L^g is the heating by 
Hen photoionization, Lh is the heating by Hi photoionization, Lcmb is the cooling off the mi- 
crowave background, Lr is the cooling by recombination, Lff is the cooling by free-free emission, 

and La is the cooling due to line excitation and collisional ionization. We have separated the terms 
for He II and H I photoionization because He II plays a dominant role for heating, but the other 
cooling terms include both hydrogen and helium. 

To sec what the important terms for the thermal balance of the IGM are, we now evaluate 
all the heating and cooling terms at the conditions where we have measured the temperature most 
accurately: a temperature T = 20000 K at a gas density A = 1.4, at 2; = 3. It is convenient 
to define the quantities T' = 2/{3kHT)L, for each subscript corresponding to every heating and 
cooling term. We assume the model ^bh^ = 0.019, CIq = 0.3, Aq = 0.7, and Hq = 65kms~^ Mpc~^. 
We evaluate first the total cooling: the dominant term is adiabatic cooling, which is equal to —2 on 
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the right-hand-side of equation (10) if we assume a rate of expansion equal to the Hubble rate (i.e., 
a constant A). This assumption is of course not exact; every gas element expands at a different 
rate, causing a dispersion in the temperature-density relation. However, at a density A = 1.4, the 
average rate of expansion is in fact not very different from the Hubble rate. In addition, when we 
consider the evolution of the temperature at a fixed A, the dT/dA term in the total temperature 
derivative of the left-hand-side of equation (10) should partly compensate for the effect of expansion 
(canceling it exactly if 7 — 1 = 2/3). 

The next most important contribution is cooling off the microwave background, which is given 
by Tqj^^ = (8(TTa^^J^^g)/(3i^mec)(ne/n) = 0.58 (where ax is the Thompson cross section, mg 
the electron mass, Ug the electron density, and n the total particle density). Notice that this term 
becomes more important at higher redshift, growing as (1 -|- zf'l'^ . We compute the other cooling 
rates using the formulae given in Black (1981). Recombination yields = 0.22, and free- free 
emission T^j = 0.08. The atomic processes of line cooling and coUisional ionization axe completely 
negligible at this low density, for the Hi photoionization rate V ~ 10~^^ s^^ that is obtained from 
the observed abundance of quasars. The total cooling rate is therefore T' = 2.88. If this temperature 
is being kept roughly constant, as indicated by the measurement we have presented here, then the 
heating terms should approximately balance the total cooling. 

To evaluate the heating rate, we first assume ionization equilibrium; we will discuss later how 
the heating from Hell can be increased if the He 11 reionization is still in progress. The heating 
term due to ionization can then be expressed in terms of the recombination rate: Ljj =< Eh > 
{aHne){nH /n), and the analogous expression for Hen , where an and nn are the recombination 
coefficient and the number density of hydrogen, and Eh is the mean energy of the absorbed photons 
minus the ionization potential. The mean energy Eh depends on the spectrum of the ionizing 
background, and wc evaluate it as follows: assuming a background intensity per unit frequency 
Jy oc u^^ from the ionization edge at to some maximum frequency = Qm^o^ and approximating 
also the photoionization cross section as aiy) oc v^^-, we find 



< Eh >= Ih 



1 - gm^^^ [3 + 3 



(11) 



where Ih is the ionization potential. We use /? = and = 4 for both Hi and Hen , which 
adequately approximates the shape of the spectrum found in numerical calculations when the 
emitted spectrum is a quasar power-law (with /3 = 1.5), and the effect of absorption by Lyman 
limit systems is taken into account (Miralda-Escude &: Ostriker 1990; Haardt &: Madau 1996). This 
yields < Eh >= 0.43/// and < Ehs >= 0.43///e- We then obtain: = 0.67, and T^^ = 1.11. 
The total heating therefore falls short to compensate for cooling by a factor ~ 1.6. 

An obvious way to increase the heating rate is to assume that the He 11 reionization is not yet 
complete; in other words, that there are patches of low-density gas in the IGM where all the helium 
is in the form of Hen . As discussed in Miralda-Escude & Rees (1994), there are two reasons why 
the heating rate is higher during the reionization, relative to the case of photoionization equilibrium. 
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The first is that the ionization rate needs to be higher simply because every He ii ion needs to be 
ionized once during the course of the He II rcionization. The second is that all the hard photons 
will now be absorbed by a random Hell ion in the IGM, up to the frequency I'm = Qm'^o where the 
mean- free-path through the Heii IGM reaches the horizon length. For the baryon density we use 
and at z = 3, and assuming also that about 50% of all the helium is in the form of Hen in the 
diffuse IGM (and not in dense clouds having a small covering factor over the Hubble length), this 
maximum frequency is given by qm = 13, i.e., a frequency 13 times higher than the Hen ionization 
edge. The mean energy of the absorbed photons is therefore equal to the mean energy of the 
emitted photons up to this maximum frequency, without weighting them with the photoionization 
cross section: 



< EHe,r >= I He 



- 1 



(12) 



where the subscript r in < EHe,r > indicates the mean energy per absorption during reionization, 
and the emitted spectrum from sources is Ji^ oc L'~f^' . For qm = 13 and f3' = 1.5, we obtain 
< EHe,r >= 0.71///e- Since the recombination rate for Hell at the mean density is equal to 1.6 
times the Hubble rate H^^ at z = 3, and if the reionization is occurring over ~ a Hubble time 
near z ~ 3, it is reasonable to expect that the additional heating rate due to Heii reionization is 
comparable to the heating rate due to balancing recombinations of Hen . We therefore conclude 
that the heating from He ii reionization can reasonably account for the IGM temperature we have 
determined here. 



7.2. Usefulness for Measuring the Baryon Density from the Lyman Alpha Forest 



One of the applications that the development of the new theory of the Lya forest based 
on structure formation has had is to provide a measurement of the baryon density through its 
effect on the mean transmitted flux. For a fixed distribution of temperature, over-density, and 
peculiar velocities in the IGM, the Lya optical depth at any point in the spectrum is proportional 
to {p,Bh?'Y H{z)~^Tz\2-, where V = lO"-*^^ r_i2 s~^ is the photoionization rate due to the cosmic 
ionizing background. To be specific, we define the parameter 



0.3 X 0.652 



1/2 



-12 



-1/2 



(13) 



[where we have used H{z) ~ Ho{l + z)^^'^^^/'^ , which is highly accurate at the relevant redshifts and 
in a flat universe]. As discussed by Hernquist et al. (1996), Miralda-Escude et al. (1996), Ranch 

et al. (1997), Weinberg et al. (1997), and Paper I , a measurement of ujb can be translated into a 
lower bound on il^/i^ by using the contribution to the ionizing background from known quasars as 
a lower bound on r_i2. This lower limit is on the high side of the range of fib that is allowed by 
primordial nucleosynthesis: fife/i^ > 0.02 (Rauch et al. 1997; Paper I ). 
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One of the main model uncertainties in deriving the relationship between the parameter lob and 

the predicted mean transmitted flux is the mean IGM temperature. A higher temperature implies 
a lower recombination coefficient, and therefore a lower neutral hydrogen density. This implies 
that in order to reproduce a given observed mean transmitted flux, the mean density Qb needs 
to be increased further to compensate the reduced recombination coefficient. The measurements 
of the IGM temperature reported here and in Schaye et al. (2000), and Ricotti et al. (2000), all 
coincide in finding temperatures that are high compared to what is expected if the IGM is heated by 
photoionization and is in ionization equilibrium. As wc have discussed, these higher temperatures 
can probably be understood as a result of the Hell reionization. independently of its cause, the 
higher temperature implies an even higher value of a;^ than was obtained previously, which we 
can easily determine by modifying the temperature in the simulation to match the observations, 
as described earlier in §5.3. We find that, when the temperature in the simulation is increased to 
match the observed one, the derived value of lub is increased slightly to lob =(0.0336 ± 0.0020, 
0.0288 ± 0.0023, 0.0248 ± 0.0017) at z =(3.9, 3.0, 2.4), from the previous values (0.0329, 0.0274, 
0.0245) when the temperatures in the simulations arc not modified (Paper I ). The errors are derived 
from the observational errors in the determination of the mean flux decrement from Paper I . 

Changing the temperature of the simulation affects the value of not only by modifying 
the recombination coefficient, but by increasing the amount of thermal broadening, which can 
spread the absorption in saturated regions to the outskirts of absorption lines, increasing the mean 
absorption for a given lob- We find that the effect of thermal broadening is less important than 
the effect of the reduced recombination coefficient. As an example, if we replace a{T) a{T + 
3000i^) in every pixel in the simulation, the inferred ujb increases by 0.0017 (where a(T) is the 
recombination coefficient), while the replacement (Jb{T) crb{T + 3000i^) changes ujb by -0.00046 
(where crb{T) is the dispersion of the Gaussian thermal broadening). Dynamical effects caused 
by the increased pressure of the gas would probably go in the same direction as the thermal 
broadening, since they would tend to spread the gas in absorption systems over wider regions. 
However, it seems unlikely that any such dynamical effects (which can only be investigated by 
running the same simulation with different temperatures) can be more important than the thermal 
broadening effect. The examples shown in Theuns et al. (1999b) (see their Figure 6) appear to 
confirm that the dynamical effects of increased pressure are not more important than the increased 
thermal broadening when the gas temperature is raised. 

With the statistical error bars we have obtained on the T — Ag relation , we can place more 
conservative lower bounds ouojb than the ones obtained in Paper I . The lowest allowed value of a;^ 
needed to account for a given observed mean transmitted flux is obtained when is minimum and 

7 — 1 is maximum in equation (1), because that yields the minimum temperature for the low-density 
gas that determines the optical depth in unsaturated regions of the Lya spectrum. We set equal 
to the measured value minus twice the statistical error bar given in Table 5, and 7—1 = 0.6, which 
is the value valid when the IGM has been in photoionization equilibrium for a long time (Hui & 
Gnedin 1997). Any uniform heating of the IGM, such as that caused by reionization, should give 



-49- 



rise to a lower 7 — 1, although shock-heating can increase 7 — 1 above 0.6, the simulations show 

that this happens only at high enough gas densities that the Lya absorption is already saturated. 
The error in our observational determination of 7 — 1 is too large to give us a better constraint 
than 7 < 0.6 (see Table 5). 

The results of this exercise are ujb >(0.0270, 0.0192, 0.0209 ) for z =(3.9, 3.0, 2.4), at 95% 
confidence, including the error from the mean flux decrement and the temperature measurements 
(added in quadrature). Using the lower bound obtained in Ranch et al. (1997) of F— 12 > 0.7 
in the range 2 < z < 3, obtained by counting only radiation from the observed quasars, and not 
including the power-law extrapolation of the quasar luminosity function that has been observed only 
at redshifts 2; < 2, we obtain fi^/t^ > 0.017. This result is still consistent with the determinations 
of the deuterium abundance (Buries & Tytler 1998). However, if the quasar luminosity function 
extends to low luminosities with a similar power-law slope as observed at z < 2, or if emission from 
galaxies increases significantly the intensity of the ionizing background, then the higher baxyon 
density implied would come into conflict with the primordial nucleosynthesis predictions and the 
observed deuterium abundance. 

In summary, we have reached the following conclusions: 

1. The temperature of the IGM is ~ 20000 ± 2000K at density 1.4 times the mean, independent of 

redshift, although an increase of ~ 3500K from z = 3.9 to 2; = 3.0 cannot be ruled out. 

2. The high temperature cannot be explained by heating in ionization equilibrium, and probably 
indicates on-going He 11 reionization. 

3. The contribution of temperature uncertainty to the uncertainty in the baryon density required 
by the observed mean flux decrement in the Lya forest is now well constrained. 

We thank Adam Steed and David Weinberg for helpful comments on the manuscript. 

A. THE PROFILE FITTER 

This algorithm has three input parameters that control how the fitting proceeds: sets the 
amount that the flux must increase from the center point to the edges of the window before a fit 
will be attempted, Wmin sets the minimum size of the window within which a fit is performed, and 
Pq sets the quality of fit that will be accepted. In this paper we set Wmin = 2. 

Three more input parameters effect the speed of the code but are not important to the results: 

Eg controls the degree of symmetry around a central pixel that is required for a fit to be attempted, 
Ec sets the level of flux decrease, from the center pixel to the window edges, at which a point will 
be eliminated from consideration for fitting, and Wmax sets the maximum allowed window size. 
These parameters are set to values large enough that they do not actually eliminate any profiles 
that would otherwise be accepted. 

Before we describe the algorithm in detail, a few more terms must be introduced: We are going 
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to fit pieces of tfie spectrum that have center point P and extend iiW pixels to either side of P. 
The width of the fitting window, W, will be adjustable but constrained to Wmin ^ W < Wmax- 
The transmitted flux at a point P is F{P). The error in the sum or difference of the flux at two 
points Pi and P2 is a{Pi,P2) = [a{Pi)'^ + a{P2)'^Y^'^ ■, where (y{P) is the observational error in the 
flux at point P. The minimum acceptable probability for is -Po (we need to define Pq by the 
probability because there will be varying numbers of degrees of freedom in the fits). 

For a given spectrum the algorithm that we use is the following (the reader should keep in mind 
that, except for the added complication of setting the window position and width, this procedure 
just fits a single Voigt profile to each absorption maximum by minimization): 

1. Scan along the spectrum pixel by pixel searching for places where \F{P—Wmin)—F{P+Wmin)\ < 
Es a{P + Wmin,P-Wmin)- Also require that F{P)-F{P±W^in) < (j{P,P±Wrrun), where the 
flux at the ±Wmin points is averaged. These places are candidates for a symmetric, non-concave 
profile. 

2. If there is a significant increase in fiux at the edges of the window, so that F{P) — F{P±Wmin) > 
Ed (j{P,P ± Wmin), go ahead and fit Equation (4) to the absorption. If there isn't a significant 

increase try to expand the window. 

3. To expand the window require that symmetry is maintained when W is increased, i.e., \F(P — 
W) — F{P + W)\ < Es a{P + W,P — W). If the window can be expanded return to step 2 to check 
if a fit can be done with the enlarged window, i.e., if F{P) — F{P ± W) > E^ a{P, P ± W). 

4. If the region can't be fit, but also can't be expanded, eliminate the candidate point. Also elimi- 
nate the point if the window size has been increased to Wmax without meeting the requirement for 
fitting. 

5. Set initial parameters for the fit using F{P) to set Tc and [F(P + W) + F{P — W)]/2 to set aj,. 
Set Vc = 0. If F{P) < set Tc = 10. 

6. Minimize using the flux values and their error bars in the range of points between P + W and 
P — W. Require that \vc\ < 0.5 pixels (outside this range is covered by other candidate points). 

7. Eliminate candidate if P{> x^, i^) < Pq- 

8. If P falls within W of a previously accepted candidate, eliminate the candidate with a smaller 
value of P(> x^)^)- This does not eliminate any independent profiles because a single Gaussian 
would not fit if the window contained multiple lines. 
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Table 1. Power-law fits to the T — Ap relation . 



z To 7-1 (for T) Tq 7 - 1 (for T) 





(K) 




(K) 




4 


14024 


0.21 


13873 


0.22 


3 


15843 


0.30 


15584 


0.30 


2 


12764 


0.57 


12954 


0.54 



Note. — r is the optical-depth-weighted average temperature at points in spectra (fitted vs. the 
density similarly averaged). 
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Table 2. 



betting 


A T 


A T 


Q 


— 95 


1 Q8 




1 3 8 






11 


1 4 3 


— 17 


248 


19 


14.1 


£^0 = 15 


264 


22 


14.3 


E^ = U 


272 


28 


14.1 


Ec = l3 


271 


34 


13.6 


Ed = 12 


281 


40 


13.5 


Ed = ll 


285 


56 


12.4 


Ed = 8 


364 


95 


12.6 


Ed = 7 


386 


137 


10.9 


Ed = 5 


474 


263 


7.8 


Ed = 12 (0.001) 


284 


41 


13.5 


Ed = 12 (0.1) 


243 


35 


12.5 


Ed = 12 (nn) 


236 


1 


15.3 


Ed = 12 (nc) 


262 


38 


12.9 


Ed = 12 (nr) 


286 


39 


13.7 


Ed = 12 (nr) 


306 


30 


15.1 


E,i = 12 (ur) 


308 


35 


1-1.7 



Note. — The quality measure Q = {Ng- Nh){Ng + iVb)"^/^, where Ng is the number of fitted 
features satisfying K < B - T < 3000 K and Nb is the number satisfying -3000 K < B - T < 
K. Entries labeled (0.1) and (0.001) have P{> x^) >0-l and 0.001, respectively [the rest have 
P{> X^) > 0.01]. The label (nn) means no noise, (nc) means no continuum fitting approximation, 
and (nr) means a new set of random numbers was used for the added noise in each example. 
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Table 3. Basic statistics of the observational data in each redshift bin. 





^max 


z 


F 


n 


pixels 


path length 


Ed 


absorption 




Tmax 


lines 














(kms-^) 




lines 






used 


3.39 


4.43 


3.89 


0.48 


0.029 


35120 


70893 


8 


281 


3.8 


47 


98 


2.67 


3.39 


2.99 


0.68 


0.011 


35283 


87308 


12 


284 


1.0 


19 


160 


2.09 


2.67 


2.41 


0.81 


0.028 


36150 


104581 


9 


223 


0.41 


5.4 


179 



Note. — The minimum (maximum) optical depth of fitted lines used in the temperature mea- 
surement is given by r^m {Tmax)- The mean noise level in the spectra is n. 
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Table 4. Binned Results of Temperature from Observations 



'^c,med 


Tc,min 


Be 


AT 


A. 






(K) 


(K) 






z = 3.9 


0.72 


0.19 


28200 ± 16700 


16900 


0.48 


1.53 


1.25 


17400 ±11200 


11100 


0.67 


3.22 


2.34 


21200 ± 4100 


6100 


0.97 


7.83 


5.97 


22400 lb 7900 


3100 


1.66 


z = 3.0 


0.43 


0.18 


25800 ± 2200 


13800 


0.61 


1.03 


0.62 


24600 ± 3700 


4400 


0.90 


1.55 


1.29 


21800 ± 2100 


3000 


1.15 


4.80 


2.93 


25600 ± 5600 


1600 


2.53 


z = 2A 


0.52 


0.18 


22600 ± 2600 


2300 


0.99 


0.96 


0.74 


23200 ± 8500 


3000 


1.44 


1.59 


1.48 


33000 ± 5400 


1000 


2.07 


3.88 


3.18 


41600 ± 12100 


-700 


4.01 



Note. — The temperature cutoff Be, systematic offset AT, and gas density Ag, are given for 
each bin with minimum optical depth Tc,min and median optical depth Teamed- 
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Table 5. Power-law fits to Temperature from Observations. 



^obs 








7-1 


A. 


n 








(K) 


(K) 






(K) 




3.89 


4 


20200 ± 2700 


20100 ± 2800 


0.43 ± 0.45 


1.42 ±0.08 


17400 ± 3900 


6.54 


2.99 


4 


20200 ± 1300 


20300 ± 1400 


0.29 ±0.30 


1.37 ±0.11 


18400 ± 2100 


1.83 


2.99 


3 


19600 ± 1500 


20400 ± 1800 


0.33 ± 0.28 


1.24 ±0.10 


18300 ± 1800 


1.74 


2.41 


3 


22600 ± 1900 


20700 ± 1900 


0.52 ±0.14 


1.66 ±0.11 


17400 ± 1900 


1.07 


2.11 


2 


23 100 ± 2000 


22800 ± 2100 


0.51 ±0.14 


1.17 ±0.10 


19200 ± 2000 


0.98 



